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We review the properties of nuclear halos and nuclear skins in drip line nuclei in the framework 
of the spherical Hartree-Fock-Bogoliubov theory with continuum effects and projection on good 
particle number with the Gogny force. We first establish the position of the un-projected HFB drip 
lines for the two most employed parametrizations of the Gogny force and show that the use of finite- 
range interactions leads almost always to small-sized halos, even in the least bound nuclei, which is 
in agreement with most mean-field predictions. We also discuss the size of the neutron skin at the 
drip line and its relation to neutron asymmetry. The impact of particle-number projection and its 
conceptual consequences near the drip line are analyzed in detail. In particular, we discuss the role 
of the chemical potential in a projected theory and the criteria required to define the drip line. We 
show that including particle number projection can shift the latter, in particular near closed shells. 
We notice that, as a result, the size of the halo can be increased due to larger pairing correlations. 
However, combining the most realistic pairing interaction, a proper treatment of the continuum and 
particle number projection does not permit to reproduce the very large halos observed in very light 
nuclei. 
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Neutron-rich nuclei present us with unique opportuni- 
ties to test nuclear models. As the asymmetry between 
the number of neutrons and protons in atomic nuclei in- 
creases, a number of new phenomena appear such as neu- 
tron skins, nuclear halos or shell-melting li!, 2]. With cur- 
rent and on-going developments of radioactive ion beam 
facilities, new territories of the nuclear chart are open 
to exploration and data from very neutron-rich nuclei is 
likely to upset the existing theories [3|] . 

The so-called nuclear skin in neutron-rich nuclei is 
caused to a large extent by the isospin asymmetry. In 
the mean-field picture of the atomic nucleus, it should 
therefore depend mostly on the iso-vector component of 
the effective hamiltonian or Lagrangian. By contrast, the 
understanding of nuclear halos in the context of mean- 
field theories is quite fragmentary and it is not clear if one 
can relate this phenomenon to a particular term of the 
effective interaction or density functional. The extreme 
difficulties to approach the drip lines in heavy nuclei, 
where mean-field theories are most often employed, also 
forbids to test the calculations against experiment. In 
light nuclei, where experimental data is available, few- 
body models that introduce a core surrounded by one 
(two-body models) or two particles (three-body models) 
are very successful [1, However mean- field models 
are notoriously unreliable in these extremely light sys- 
tems unless severe corrections beyond the mean-field are 
included. 

In heavy nuclei it is commonly thought that nuclear ha- 
los should be interpreted as resulting from the coupling 
to the continuum via residual interactions such as pair- 
ing correlations [GJ . The spatial delocalization of contin- 
uum states gives a simple motivation for such interpreta- 
tions. Nevertheless mean-field theories show significant 



variations in their predictions of halos. Several major 
difhcultics can explain these discrepancies: firstly, the 
neutron drip line is not known beyond the Oxygen el- 
ement. This lack of experimental data is going to be 
partially filled in the near future but imply that models 
can not really be benchmarked against experiment. Sec- 
ondly all self-consistent approaches to nuclear structure 
rely on the parametrization of some effective interaction, 
Lagrangian, or energy density functional. The extrapola- 
bility of such interactions to regions of very large neutron 
excess is by no means guaranteed by the theory. Fur- 
thermore, beyond mean-field effects such as symmetry 
restoration mechanisms or configuration mixing might be 
playing a different role in these extreme regions as in the 
valley of stability. 

In a recent article, we proposed a simple and ef- 
fective method to treat the continuum in configuration 
space for mean-field theories based on finite-range inter- 
actions of the Gogny type. We showed that our procedure 
provides the correct asymptotic behavior of the nuclear 
wave-functions and established the proton and neutron 
drip lines with the DIS interaction. Restoration of the 
particle number symmetry was also discussed and it was 
found that in some cases, in particular when the neutron 
number is near a closed shell, the variation after projec- 
tion (VAP) method could shift the position of the drip 
lines by 2 neutrons. In this paper, we apply our method 
to the specific problems of nuclear skins and nuclear ha- 
los near the neutron drip line. We will furthermore re- 
view several cases that have been proposed in the past, 
either in the framework of the Skyrme Hartree-Fock Bo- 
goliubov theory or the relativistic mean-field and discuss 
the impact both of finite-range interaction and symme- 
try restoration effects. To our knowledge this constitutes 
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the first attempt to address this problem in explicitly 
symmetry-conserving mean-field approach. In section I 
we briefly recall the main features of our approach as 
well as the Helm model that has been traditionally used 
to describe nuclear halos. In section II we discuss the 
nuclear skins and halos in the standard framework of the 
HFB theory with the Gogny interaction. In section III 
we investigate the impact of particle-number symmetry 
restoration on halos. 



I. BRIEF DESCRIPTION OF THE METHOD 

In the nuclear mean-field approach, the energy of the 
nucleus is computed as the expectation value of a two- 
body Hamiltonian on a trial wave- function 8]. Besides 
the relativistic approaches [!|, there exists two main fam- 
ilies of two-body effective forces to this date, the zero- 
range Skyrme Il0l| interaction and the finite-range Gogny 
one [11]. Both are empirical effective forces and there ex- 
ists a number of realistic parametrizations. Skyrme forces 
lead to a local energy density which is the basic building 
block of the nuclear Energy Density Functional (EDF) 
theory The Gogny interaction, because of its finite 
range, is non-local and computationally more involved, 
for this reason the calculations are most conveniently car- 
ried out in configuration space, i.e. the solutions to the 
Hartree-Fock (HF) or Hartree-Fock-Bogohubov (HFB) 
problem are expanded on a given basis. 

The harmonic oscillator (HO) basis has always played 
a special role in configuration space calculations, as its 
eigenfunctions are given analytically and are separable. 
However, a well-known deficiency of this basis is that it 
is made exclusively of bound-states since the underlying 
potential has infinite walls. In practice, since calcula- 
tions are always performed in a given truncation scheme 
(for a fixed cut-off of the basis) the localization of all 
HO basis states imply that the physical wave-functions 
of the system will always acquire a Gaussian asymptotic, 
including the weakly-boimd and positive-energy states. 
This is clearly unrealistic, as spherical continuum states 
should be spherical waves. The consequences of this de- 
ficiency become more serious close to the drip line, as 
pairing correlations can couple discrete bound-states to 
the continuum. It is thus critical to properly describe 
the continuous spectrum even at the level of ground-state 
calculations p^ . 

There exists a number of techniques to take into ac- 
count the continuum in nuclear structure calculations, 
and it is not the purpose of this article to discuss in de- 
tail the merits of every one of them. Let us just men- 
tion briefly for completeness: coordinate-space Skyrme 
HFB theory with either vanishing 14] or outgoing-wave 



boundary conditions [15|, coordinate-space Relativistic 
Hartree-Bogoliubov with finite-element method [l^ , the 
use of the Gamow basis in the Skyrme HFB theory [13] , 
in the Continuum Shell Model [18| and in the Gamow 
Shell Model (ISJ. At the present time technical difficul- 



ties in including the full continuum with the exact reso- 
nant and non-resonant spectrum lead to the consequence 
that the most advanced theories are only applied with 
simple model interactions that are tailored to capture 
the main physical properties of the system. Only in the 
coordinate-space HFB approach realistic Skyrme interac- 
tions were employed with density-dependence zero-range 
forces in the pairing channel (requiring the introduction 
of a cut-off in the quasi-particle spectrum or a regulariza- 
tion procedure [lO])- Moreover, as far as mean-field based 
theories are concerned, no attempt has been done to in- 
clude with the coupling to the continuum the restoration 
of broken symmetries or collective motion. 

Therefore, in order to combine the flexibility of con- 
figuration space calculations with the necessary inclu- 
sion of the continuum, it has been proposed in 0, [2l| 
to work in a basis made of the eigenstates of the Woods- 
Saxon potential. The latter are obtained by integrating 
the Schrodinger equation in a box of size Rtox with a 
mesh size h. In practice Rtox — 20 fm and h — 0.1 fm 
are sufficient to obtain a good convergence of the solu- 
tions. Boundary conditions are set on the walls of the 
box. As usual, several choices are possible. Outgoing 
wave boundary conditions lead to wave- functions that 
are not square-integrable, and special techiiiques must 
be employed to overcome this difficulty [1^ [Ij, [13) ■ 
Vanishing boundary conditions guarantee that the basis 
functions are square-integrable and can thus be normal- 
ized at the price of eliminating all the continuum states 
that do not have a node on the walls of the box. It 
was shown in [2^ that both techniques essentially lead 
to very similar results as far as bound-states and bulk 
properties of nuclei are concerned. In the following, we 
use vanishing box boundary conditions. 

In our calculations we use the finite-range Gogny inter- 
action [ll| . The same interaction is used in the particle- 
hole channel (mean-field) and particle-particle channel 
(pairing) and both the direct and exchange contributions 
coming from all the terms of the interaction are taken 
into account in the calculation. The finite-range of the 
force in the pairing channel allows to avoid the divergence 
problem (in momentum space) and cut-off dependence of 
zero-range forces. All of our calculations are performed 
in spherical symmetry. 

To obtain quantitative information on the neutron halo 
in neutron-rich nuclei we make use of the Helm method 
[13, [H,!!!, 1131. Firstly, the neutrons (protons) form fac- 
tor is computed as the Fourier transform of the neutrons 
(protons) density. In spherical symmetry, this leads to: 
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where q is the momentum, jo is the spherical Bessel func- 
tion of order and Prir) is the density (r standing for 
neutron or proton). This form factor built out of the re- 
alistic one-body density, in our case calculated with the 
Gogny force, is then compared to the Helm form-factor 
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obtained from the convolution of the Gaussian profile 
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with a sharp density profile: p{r) — po for r < Rq and 
p{r) = elsewhere. Since this model is presented in 
details in the references quoted, we simply recall the for- 
mulas we are going to use. The two parameters Rq and 
a of the model are determined in the following way. 

The rms radius Rrms is defined as the squared root of 
the mean- value of the operator f^. It is extracted from 
the nucleonic density: 



Rr 



\ 



p{r) 
cPrp{r) 



(3) 



For the Helm radius one straightforwardly obtains 
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where Rq is the diffraction radius: 

Rq = 4.49341/gi 



(4) 
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and qi is the first zero of the realistic form-factor ([T]) ob- 
tained in our theoretical approach. The surface thickness 
a is defined as: 
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where N is the number of particles, ji the spherical Bessel 
function of order 1 and q,n is the first maximum of the 
realistic form- factor ([1]). 

At this point we should note that the method does 
not carry out any information on the eventual decorre- 
lation between a core and a few valence particles. It 
only provides a simple and fast method to assess the 
spatial extension of the nucleus and an excellent starting 
point to determine the best halo candidates. However, in 
few-body nuclear models, the nuclear halo is often inter- 
preted as one single nucleon or a pair of nucleons orbit- 
ing around a core, see e.g. [3l|, [32| for two-body models 
and [H, m, |3^ for 3-body models. In order to reconcile 
these cluster approaches with a mean-field description 
of the nucleus, a more detailed analysis of the density 
should be carried out. Alternative techniques have been 
proposed to cure this deficiency [3(|. 

It is usually convenient to multiply the rms and the 
Helm radius by •y/5/3. The quantities 
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are related to the underlying shape of the nucleus. A 
measure of the nuclear halo is then provided by the quan- 
tity: 



The neutron skin can be defined in various ways depend- 
ing on which type of radius is considered. Its general 
expression is: 
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where R can be either the geometrical radius, the Helm 
radius or the diffraction radius. It was argued in [s^] that 
the best approximation to the neutron skin is obtained 
when taking the Helm radius, as the latter is somewhat 
rid of spurious contributions coming from the neutron 
halo. 



II. NUCLEAR SKINS AND HALOS WITH 
FINITE-RANGE INTERACTIONS 

Systematic calculations near the neutron drip line have 
been carried out using the spherical HFB code in the 
Woods-Saxon basis that was presented in [7]. The basis 
was constructed from the eigenstates of the WS potential 
with the universal parametrization of [13] applied to the 
Z = 126 and N — 184 nucleus. The Schrodinger equation 
was integrated in a box of Rbox = 20 fm with vanishing 
boundary conditions. All eigenstates with i < 15 and 
n < 18 were retained in the basis. As was shown in 
7], such a choice guarantees a good convergence of the 
subsequent HFB calculation. 



A. Determination and properties of the neutron 
drip line 

There exist few parametrizations of the Gogny interac- 
tion: in our calculations wc considered the parametriza- 
tions Dl of [m and DIS of [38|. For each of them the neu- 
tron drip line was calculated based on the requirement 
that the one-neutron Sn = B{N, Z) - B{N - 1, Z) sep- 
aration energy must be negative for bound nuclei. Since 
S2n = Sn + Sn-1, the Criterion Sn < is stricter than 
the condition that the two-neutron separation energy is 
negative. In the HFB theory the one-neutron separation 
energy Sn is approximated by the neutrons Fermi en- 
ergy A„ = dE/dN Ki —Sn- A nearly equivalent condition 
to define the one neutron drip line is therefore A„ > 0. 
When HFB pairing correlations vanish (case of closed 
shells) , the value of the chemical potential A is meaning- 
less and can not be used to define the drip line any more 
(Hartree-Fock limit). However, in the HF approach and 
within the approximation of the validity of Koopmans' 
theorem [39], the stability of a nucleus is simply gov- 
erned by the position of the last occupied level: if it has 
positive energy, then the nucleus is unbound with respect 
to particle emission. 

We display in Table U the one neutron drip line nu- 
clei obtained with both interactions. In the presence of 
pairing correlations the criterion A„ > has been used. 
For the neutron shell closures = 82 (DIS: elements 
Z = 36 to Z = 40, Dl: elements Z = 36 and Z = 38), 
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N = 126 (DIS: elements Z = 52 to Z = 64, Dl: ele- 
ments Z = 54 to Z = 62) and N = 184 (DIS: elements 
Z = 80 to Z = 92, Dl: elements Z = 80 to Z = 90) 
the neutron pairing correlations vanishes and we have to 
rely on Koopmans' theorem. The columns correspond- 
ing to the DIS interaction were already presented in Q 
and are recalled for comparison. We would like to com- 
ment at this point that, due to some technical problems 
with our previous codes at the aforementioned shell clo- 
sures, in [3] the drip line was predicted with two neu- 
tron less for the elements Kr, Te, Xe, Ba, Hg and Pb. 
We also show in Table |T] the difference in the number of 
neutrons between the drip line nuclei with the Dl and 
DIS interaction: AiV = NdripUne{Dl) ~ NdripUne{DlS). 
In general the Dl parametrization predicts a drip line 
with more neutrons, probably due to the fact that it pro- 
vides more pairing correlations than the DIS one. Let us 
emphasize that all calculations performed in this work 
are restricted to spherical symmetry. Several of the nu- 
clei listed in Table U may be deformed in their ground- 
state 4l|. Symmetry- unrestricted HFB calculations of 
neutron-rich nuclei would most likely shift the position 
of the drip line in several places. 



TABLE I: Table of spherical HFB one neutron drip line nuclei 
obtained with the DIS and Dl interactions. The columns 
marked AA'^ represent the shift of the drip line (in number 
of neutrons) when using the Dl interaction compared to the 
DIS. 
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For each interaction, the quantity SRhaio of Eq. ([5]) 
was then computed at the drip line, i.e. for each ele- 
ment listed in Table [J The results are shown in Fig. [T] 
We find a downward trend of SRhaio as a function of Z 
superimposed with oscillations. Both features are well 



understood. The decreasing behavior has to do with 
the well-known fact that light nuclei have larger halos. 
The oscillations are related to the neutron magic num- 
bers: To the five minima (for the DIS, for example, 
Zmin — 10, 22, 40, 64 and 92 ) correspond to neutron num- 
ber 20, 50, 82, 126 and 184, see TableU For proton num- 
bers two (four or six) units larger than a given Zmin a 
few neutrons occupy a new large j-shell thereby inducing 
pairing correlations and producing a halo. For Z values 
much larger than a given Z„iin the number of neutrons in 
the shell becomes large and the halo disappears. 




40 60 80 100 

Proton Number Z 

FIG. 1: Measure of the neutron halo: SRhaio ~ Rgeom — 
Rneim in fm for spherical Gogny HFB calculations in the 
WS basis with the DIS (plain squares) and Dl (open circles) 
interactions. 

As noticed in [s^l, the size of the halos is correlated 
with the corresponding chemical potential of the HFB 
solutions: Larger halos correspond to nuclei with values 
of A„ close to zero and smaller ones to large A„ values. As 
far as the effect of the parametrization of the Gogny force 
is concerned, we observe that the largest difference takes 
place in ^^Mg which is not bound with DIS interaction 
while it is bound for the Dl interaction. Apart from this 
particular nucleus, both parametrizations of the Gogny 
force give very similar results, even though the isotopes 
of the drip line elements are sometimes very different, for 
example ^-'^^Er with the Dl interaction and ^oeg,. ^fJ{l\^ 
the DIS. 

It is instructive to compare our results with the work 
of [lOl- It was pointed out in this reference that the 
size of the halo, as measured by the quantity SRhaio, 
significantly depends on the interaction used. A similar 
conclusion was reached in [36l | using a slightly different 
analysis procedure. In our case both parametrizations 
provide rather similar results in spite of the fact that the 
numerical values of the DIS and Dl parametrizations are 
quite different. It is also interesting to note that both 
parametrizations can lead in some cases to significantly 
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different drip lines - in the case of the Dl interaction for 
example, the drip line near Palladium isotopes (Z = 46) 
and Erbium {Z = 68) and Ytterbium [Z = 70) is located 
8 and 10 neutrons further away, respectively, than in the 
case of the DIS interaction. Yet, as mentioned, the size 
of the halo remains very similar. 
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FIG. 2: Upper panel: Neutron Rgsomin) and RHeimin) and 
proton Rgeomijp) and RHeim{p) radii for the Ni isotopes calcu- 
lated with the DIS Gogny interaction. Lower panel: Neutron 
and proton halo parameter SRhaio along this isotopic chain. 

Figure [1] shows which elements can be considered as 
the best halo candidates. For each such candidate, the 
inspection of the full isotopic sequence from drip line to 
drip line can provide information on the swelling of the 
nuclear skin and the transition skin to halo. As a first ex- 
ample, we show in Fig. [2]the case of the isotopic chain for 
Nickel element. For the DIS interaction, this element has 
one of the largest halo. Furthermore, the same isotopic 
line was also studied in the framework of the Skyrme- 
HFB (SLy4 and SKP interactions) and RHB (NLSH and 
NL3 lagrangians) theories [SC^, which therefore gives us 



results from three different sorts of mean-fields. In the 
upper panel of figure [2] both the geometrical and Helm 
radii are plotted for the neutron and proton along the Ni 
isotopic chain. In the lower panel, the halo parameter 
SRhaio is plotted for the neutrons and protons. In nu- 
clei far away from the drip line, the difference between 
geometrical and Helm radius is very small, reflecting the 
neglegible coupling to the continuum near the valley of 
stability. At iV = 50 we observe the last shell closure 
and immediately after the onset of pairing correlations, 
which translates into a rapid increase of the halo param- 
eter. As was pointed out in [3(| , a shortcoming of the 
Helm method is that in some cases, the quantity SRhaio 
is non-zero even in the middle of the valley of stability. 
This appears clearly in Fig. [2] for the protons along the 
entire isotopic chain. 
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FIG. 3: Upper panel: Neutron Rgeom{n) and RHeimin) and 
proton Rgeomip) and Rneimip) radii for the Sn isotopes calcu- 
lated with the DIS Gogny interaction. Lower panel: Neutron 
and proton halo parameter SRhaio along this isotopic chain. 

Another important remark is that both parametriza- 
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tions of the Gogny force tend to give "compact" nuclei, 
with relatively small halos in agreement with those ob- 
tained in [s^l with the Skyrme SkP and relativistic mean- 
field NLSH and NL3 parametrizations and in contrast to 
the Skyrme/SLy4 interaction. The case of Tin isotopes is 
even more enlightening. For this particular element, the 
position of the drip line is nearly identical in spherical 
HFB calculations with Gogny/DlS and Skyrnie/SLy4 in- 
teractions, which facilitates the comparison. In the upper 
panel of Fig. [3] we plot the neutron and proton geometri- 
cal and Helm radius from the proton to the neutron drip 
line. Note that near the neutron drip line, the halo is 
only about 0.15 fm while Skryme/SLy4 results reported 
in [30| indicate a size of about 0.8 fm. 

It could also be tempting to apply our method in some 
of the experimental cases of nuclear halos. However, as 
hinted in the introduction, we are faced with one major 
difficulty: most of the halo candidates are very light nu- 
clei with Z < 6 at the drip line like ^^Li, ^''Be and ^^B. 
For all the elements with Z < 8 the experimental drip 
line is rigorously established, in the sense that isotopes 
bey ond the drip line are proved to be particle-unstable 
pa | . The application of our spherical Gogny-HFB calcu- 
lations, whether the Dl or DIS interaction is used, gives 
the correct drip line isotope for Lithium (see next section) 
but fails to reproduce the experimental data for elements 
B, C and O. Three main mechanisms, possibly combined, 
can be the source of this discrepancy: (i) the interactions 
used are not extrapolable in these light nuclei (ii) addi- 
tional mean-field symmetries must be broken, e.g. rota- 
tional invariance (iii) correlations beyond the HFB level 
must be included. It is almost certain that the fit of the 
interactions can be improved, but it is today difficult to 
assess to which extent this would affect the predictions of 
nuclear halos in very light nuclei. Similarly, it is not very 
clear at the moment how halos are formed in deformed 
nuclei. 



B. Discussion on giant halos 

Since our description contains the main ingredients for 
a proper description of the halo phenomenon, namely, 
a good pairing force, the incorporation of the contin- 
uum and eventually particle number projection (in the 
VAP approach) indispensable in a weak pairing regime 
we can confront our model with recent spectacular pre- 
dictions about the existence of several giant halos in light 
and medium-mass nuclei. In Neon isotopes, spherical 
coordinate-space Relativistic Hartree-Bogoliubov (RHB) 
calculations predicted that giant halos could develop for 
a number of neutrons around 30 [4^. In our spherical 
Gogny-HFB calculations the drip line is positioned at 
= 20 for both Dl and DIS interactions, which falls 
a bit short of the last known bound Neon isotope at 
iV = 24 and references therein. When deformation 
is included in the calculation, the position of the (cur- 
rent) drip line shifts to iV = 24 [41]. The very stretched 



drip line reported in (40| is somewhat surprising since 
the pairing channel was treated by using the DIS finite- 
range interaction. Moreover, up to = 20, RHB results 
for the r.m.s. radii are very similar to ours: for example 
at A^ = 20 we find a neutron r.m.s radius of r„ = 3.39 fm 
with the DIS interaction, to compare with RHB results 
of r„ w 3.42 fm. 
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FIG. 4: (color on-line) Upper panel: Two-neutron separa- 
tion energy for the Zr isotopic chain with the Dl and DIS 
parametrization. Lower panel: Neutron single-particle levels 
in the canonical basis for Zr isotopes (DIS interaction). The 
bullets represent the position of the Fermi level. 

The application of the spherical coordinate-space 
RHB, with a zero-range density-dependent force in the 
particle-particle channel, led to another prediction of gi- 
ant halos in Zirconium isotopes [4^ . Similarly as in the 
case of Neon isotopes, such predictions are rooted in the 
existence of a very stretched drip line at A^ = 100 corre- 
sponding to element ^""^Zr. Results showed in TableUand 
upper panel of Fig. |4] show that the drip line in spher- 
ical Gogny-HFB calculations is at A^ = 82 for the DIS 
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FIG. 5: Neutron Rgeomin) and RHeim{n) radii for the Zr iso- 
topes calculated with the DIS (plain symbols) and Dl (open 
symbols) Gogny interaction. 



and TV = 84 for the Dl interaction. Deformed Gogny- 
HFB calculations with the DIS interaction also suggest 
a drip line at = 82 [JJ- These results are in agree- 
ment e.g. with Skyrme HFB calculations with the SLy4 
interaction [i^ which predict the drip line at = 84. 
Other parametrizations of the Skyrme interaction have 
slightly more extended drip lines, at N = 92 for SKP 
and N = 94 for SKM* For the Gogny interaction, 

isotopes with > 82 (DIS) and iV > 84 (Dl) arc un- 
bound with respect to two neutron emission, see upper 
panel of Fig. [H Beyond drip line HFB calculations, al- 
though not realistic, can be pedagogical: in Fig. [5] we 
show the evolution of the neutron geometrical and Helm 
radius beyond the drip line for both the Dl and DIS 
interactions. As we increase the number of neutrons, de- 
localized orbitals corresponding to discretized continuum 
states become occupied and cause a very marked increase 
of the neutron radius. 

If we restrict ourselves to physical solutions at the HFB 
level we find very small halos: SRhaio ~ 0.06 fm for the 
DIS interaction in ^^^Zr and SRhaio ~ 0.11 fm for the Dl 
interaction in ^^^Zr. One may be tempted to attribute 
this small value to the collapse of pairing correlations 
which occurs at = 82. This collapse of pairing cor- 
relation can be inferred from the lower panel of Fig. |4l 
where the gap between the (occupied) hii/2 and (empty) 
2f7/2 orbital is very large. However, particle- number pro- 
jection before variation applied to this nucleus provides 
the same drip line Zr isotope and leads essentially to the 
same value of SRhaio even though pairing correlations do 
not vanish any more. 

It should be noted that our results agree with previous 
works from [H (RHB), [H] (SLY4, SKM* and SKP) as 
far as the main features of the shell structure of Zr iso- 
topes are concerned, cf. for example the neutron single- 
particle levels in the canonical basis. Fig. 1 in [4^] and 



Fig. [4] in the present work. In particular, the inflection 
point in the neutron radius at A" = 82 is reproduced 
by all models. However, all three realizations of the nu- 
clear mean-field differ as to the exact location of the drip 
line for Zirconium element, and it is this uncertainty that 
causes the widely different predictions of halo sizes in this 
particular element. 
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FIG. 6: Upper panel: Neutron skins for the Ni isotopes cal- 
culated with the DIS Gogny interaction and expressed as 
the difference of geometrical radii (plain squares) , Helm radii 
(open circles) and diffraction radii (plain triangles). Lower 
Panel: Same figure as the upper panel for Sn element. 



C. Neutron skins 

One of the main features of neutron-rich nuclei is the 
development of neutron skins as the asymmetry between 
the number of neutrons and protons increases. The 
method that we developed to include the continuum in 
our calculations allow us to compute neutron skins up to 
the drip line. As an illustration, we display in Fig. ([6|) 
the neutron skins calculated from the geometrical. Helm 
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and diffraction radius for the two isotopic chains of Nickel 
and Tin. 

As expected all three definitions of the neutron skin 
give a smooth increase with the neutron number. As 
noticed in (s^l, however, the skin calculated from geo- 
metrical radii shows a clear inflection point at iV = 50 
(Ni isotopes) and iV = 82 (Sn isotopes), which is directly 
related to the one marking the appearance of the neutron 
halo, cf. Fig. [3l By contrast, the neutron skin calculated 
from the Helm radius is a more regular function of the 
neutron number. Interestingly, although the size of the 
halo with our Gogny/DlS interaction is markedly smaller 
than with e.g. Skyrme/SLy4, the values for the neutron 
skin are much closer: In ^^°Sn, /^.R^jeim ~ 0.57 fm for 
DIS and ARHeim ~ 0.70 fm for Skyrme/SLy4 (a similar 
number is also obtained in Skyrme/SKP), cf. \30i] . 

It is instructive to compute the neutron skin for all the 
elements located at the drip line. In the upper panel of 
Fig.® we plot the neutron skin for the nuclei listed in 
Table |T1 We find an oscillatory behavior relatively similar 
to the one found for the quantity SRhaio plotted in Fig. 
[TJ In both cases, halos and skins, these oscillations can 
be somewhat related to neutron magic numbers, but the 
underlying physics is quite different. 

Skins are defined as the difference between the neutron 
and proton radius. Therefore, variations in the shape 
of the skins measures the relative increase or decrease 
of neutrons versus protons. At a neutron shell closure, 
one can add several protons without changing the po- 
sition of the neutron drip line, i.e., as a function of Z, 
the proton radius increase and the neutron one remains 
constant, which produces a decrease in the neutron skin. 
This effect is very clearly seen in Fig. [T] the minima 
at Z = 40, Z = 62 - 64 (Dl-DlS) and Z = 90 - 92 
(Dl-DlS) correspond to the last isotone with (magic) 
neutron number iV = 82, iV = 126 and N = 184 re- 
spectively. Once beyond the neutron magic number, the 
neutron radius increases very rapidly, and this translates 
into a quick increase of the neutron skin for the next few 
elements. This sharp rise is also visible in Fig. [7] in the 
range 40 < Z < 50, 62 - 64 < Z < 72. 

In fact, the oscillations of the neutron skin can be cor- 
related very neatly to the quantity (N-Z)/A. While the 
neutron excess N-Z increases with the mass number A 
(an upward trend not observed in Fig. [7]), the ratio (N- 
Z)/A fiuctuates around some average value of 0.37. In 
the lower panel of Fig. ([7]) we plot (N-Z)/ A as a function 
of Z at the drip line for the two interactions DIS and Dl 
considered in this study. We observe that the maxima 
and minima of the neutron skin correspond almost ex- 
actly to the maxima and minima of (N-Z)/A, especially 
for heavy nuclei. In light nuclei, this correspondence re- 
mains, although it is a little less obvious. We should like 
to stress that the quantity (N-Z)/ A is a direct measure 
of the ratio between the iso-vector and iso-scalar (inte- 
grated) densities. Neutron skins could therefore prove 
particularly useful to obtain experimental constraints on 
the corresponding terms of the interaction/functional. 
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FIG. 7: Upper panel: Neutron skins along the neutron drip 
line calculated with the Gogny interaction and expressed as 
the difference of Helm radii with the DIS (plain squares) and 
Dl (open circles) interactions. Lower panel: Quantity (A'' — 
Z)/A at the drip line for the DIS (plain squares) and Dl 
(open circles) interactions 



It also follows from this observation that we do not 
observe for the neutron skins the clear downward trend as 
function of Z that was observed for the halos, cf. Fig. [T] 
Neutron skins are rather a mass independent observable, 
which implies that the skin in a very light nucleus such 
as, e.g. Si {Z — 14) is of comparable size as the skin in W 
{Z — 74). The amplitude of the oscillations is also much 
smaller for the skins than for the halos, refiecting the fact 
that the (N-Z)/ A ratio does not vary too much along the 
neutron drip line. Also, the main maxima for halos and 
skins do not exactly coincide: For the DIS interaction, 
for example, the skins peak a,t Z = 14, 32, 50 and 74 and 
the halos at Z = 16, 28, 44 and 74. 
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III. INFLUENCE OF 
SYMMETRY-RESTORATION ON NUCLEAR 
HALOS 

In this section, we discuss another mechanism that can 
affect the position of the drip hne, namely the restora- 
tion of broken symmetries. We focus on the projection on 
good particle number before variation and examine sev- 
eral of its conceptual as well as practical consequences. 
The method we use to include continuum effects into our 
description of weakly-bound nuclei is indeed particularly 
suitable to include extensions beyond the mean-field. 



A. The RVAP approach 

In Q we briefly described how we can simulate the 
Variation After Projection (VAP) of the HFB solutions 
by means of the restricted- VAP (RVAP) method. Since 
along the drip lines some conceptual difhculties may arise 
we discuss the method a bit more at length. To illustrate 
how the RVAP method works we shall assume a generic 
two body Hamiltonian 



H 



Iq 



Iql'q' 



t t 



(10) 



with viqi>q> the antisymmetric matrix element 

Vlql'q' = Vlql'q> - Vlqq'l' (11) 

and (c|,Ci) the single-particle creation and annihilation 
operators in a given basis. Given the most general Bo- 
goliubov transformation 



4 



VikCi 



(12) 



the HFB method provides the product wave function 

m = E^'^\-) (13) 



that minimizes the expectation value of the Hamiltonian 
H. The matrices U and V that fix the Bogoliubov trans- 
formation of Eq. (I12p are determined by minimization of 
the functional 



e'hfb m] = 



- XnN - XzZ\<^) 



(14) 



with Ajv and Xz the Lagrange parameters that adjust the 
average number of neutron and proton. 

It can be shown ^ that the minimization of Eq. (fH|) 
amounts to the diagonalization of the matrix 



h' 

-A* 



A 

-h" 



Uk 
Vk 



-— Eh 



Uk 



(15) 



with Ek the quasi-particle energies and h' = t + T — X^ — 
Xz ■ The Hartree-Fock field F and the pairing field A are 
given by 



P;;' = y^^Vlql'q'Pq'q 
qq' 



(16) 
(17) 



with p the density matrix and n the pairing tensor defined 
by 

Pw = ($|4q|$) = 

nw = ($|q'q|$) = {V*U^)^^, . (18) 
The particle number projected energy is given by 



with P the particle number projector and 



(19) 



(20) 



To avoid cumbersome formula we do not distinguish in 
Eq. ()20|) between protons and neutrons. The simplicity 
of projection techniques lies in the fact that, while |$^) is 
a correlated many-body wave function, the intrinsic wave 
function [$) remains a product wave function, i.e. the 
variational parameters to be determined are the matrices 
U and V of Eq. ([72]) . 

In the Variation After Projection (VAP) approach the 
projected energy , see Eq. (fT9)) . is minimized directly. 
In the Projection After Variation (PAV) approach the 
HFB energy E'^p^, see Eq. ([T4|) . is minimized first and 
the projection is carried out on the HFB wave-function 
after convergence. The difference is clear: in the VAP 
method we minimize the energy of the one nucleus (Z,N) 
we are interested in while in the PAV, the energy of a 
superposition of nuclei with numbers of particle Z and N 
around the actual values. Though the variational param- 
eters are the same, the solution of the VAP equations is 
numerically much more involved than the PAV one. In a 
strong pairing regime the PAV solution might be a good 
approximation but in the general case, and in particular 
along the drip lines, the VAP one is much better. With 
finite range forces the solution of the VAP equations is 
rather involved, see [i^. Considering the additional dif- 
ficulties inherent to a proper treatment of the coupling 
to the continuum, it is clear that a full VAP solution is 
beyond the actual numerical capabilities. 

A way out of this problem is the Restricted VAP. 
In the VAP method the whole Hilbert space associated 
to the transformation Eq. (|12p is scanned in the varia- 
tional procedure. In the RVAP approach, however, only 
a restricted variational space of highly correlated wave- 
functions is allowed. In our case, since we are interested 
in pairing correlations, our restricted space should con- 
tain a whole set of paired wave-functions |$((5)) which 
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parametrically depend on the real number 5. To gen- 
erate such wave-functions with different pairing content 
and that simultaneously are consistent with our Hamil- 
tonian, we proceed in the following way: Instead of iter- 
ating Eq. [15] together with Eqs. (|16ll7p as in the usual 
HFB case, we now iterate 



h' S-A 

-S-A* -h'* 



Uk 
Vk 



El, 



Uk 
Vk 



(21) 



together with Eqs. (|16I17[) until the convergence is 
achieved. The matrices U{5) and V{5) obtained in 
this way determine the wave- functions |$((5)}. Perform- 
ing the same procedure for different 5 values we gen- 
erate the restricted correlated Hilbert space. We then 
project these wave-functions onto good-particle number 
and obtain a family of particle-number projected states 
= P^|$((5)), for 5 = l.Q,...,5max. The range 
of values for 5 is chosen in such a way that at least sev- 
eral |$^(^)) wave functions correspond to highly-paired 
states. We can then take the expectation value of the 
Hamiltonian with this set of wave- functions, i.e, using 
eq. ([T9l) . This gives us a curve [6) where, at each 
point 5, the particle number is conserved. The variational 
principle guarantees that such a curve has a minimum, 
which approaches the VAP result [i^l ■ 

To illustrate the procedure with a numerical applica- 
tion in Fig. [8] we display the un-projected energy 



E^^'^iS) 
and the projected one 



E 



PNP 



imim) 



{<P{6)\HP^\<P{5)} 



(22) 



(23) 



for the drip line nucleus ^^Ca with the DIS interaction. 
When computing the density-dependent contribution to 
the projected energy E^^^{S), the projected density 
pPNP Y^g^ been used (prescription 1 in [iq). 

Since the HFB self-consistent minimum is obtained, by 
definition, &t 5 — 1 for E^^^{5) we expect a parabolic 
behavior around this value for increasing or decreasing 
5 values. Concerning E^^^{5), at (5 = 1 projecting the 
HFB solution onto good particle number lowers the en- 
ergy and for increasing pairing correlations, i.e., values 
of 5 larger than 1, we first observe a decrease of the pro- 
jected energy up to a minimum around 5 — 1.12 followed 
by a rapid increase. Obviously the solution of the RVAP 
approach is |$((5 = 1.12)}. 

In a PNP approach the drip lines are defined in terms 
of the projected separation energies, i.e. in terms of S'j^ = 
B^iN,Z) - B^~\N - l,Z) and = B^{N,Z) - 
B^-^{N - 2, Z). In the HFB approach Sn « -A„ and 
the one neutron drip line can be easily calculated. This 
approximation is not valid any more in a projected the- 
ory and must be explicitly calculated. Since for the 
moment we are not able to project on an odd number of 
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FIG. 8: Intrinsic (plain squares) and projected (open circles) 
energy as a function of the parameter 5 for the nucleus ®^Ca. 



particles we cannot calculate the one neutron drip line 
and will therefore focus on the 2-neutron drip line. In 
the next section we discuss the meaning of A and other 
quantities in the particular context of a particle number 
projected theory. 



B. On the RVAP approach and the number of 
particles of the intrinsic wave function 

Let |<I>) be a HFB wave- function, i.e. a particle number 
symmetry violating wave function. We will now show 
that the particle number projected energy is invariant 
under transformations that change the particle number 
of the underlying HFB wave function. We define 



||.) = 6""^^!$) 



(24) 



with AN = N - No, No = ($|iV|$) and a is a real num- 
ber and we assume that (<I'|<&) = 1- The wave function 
$) can be written as [s] 



1$) 



E 

f3',N' 



c. 



1/3', TV') 



(25) 



where |/3',A^') is an eigenstate of TV with particle num- 
ber N' and [3' stands for all other necessary quantum 
numbers. The transformed wave function reads: 



1$) 



^ Cp.MP'^N') 

P',N' 

13', N' 



(26) 
(27) 
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The projected energy is given by 



E 



(l>|P^|$) 
j:f^.,P'C*^".NCp',N{P"N\H\P'N) 
Y.p»,p,C*p.>,NCp'.N{P"N\l3'N) 



(28) 
(29) 
(30) 



in an obvious way. 

The wave functions |$) and |<I>) have different number 
of particles on the average. This can be easily shown 
assuming the parameter a small enough, in this case 



($|Ar|$} 



= A^o + 2a($|(A7V)2|$). 



(31) 



up to terms. Since |<i>) is per definition a symme- 
try violating wave function, ($|(AiV)^|<I>) ^ 0, the wave 
functions |$) and |<f>) do have on the average different 
number of particles. 

We have therefore demonstrated that we can change 
the average number of particles of the intrinsic wave func- 
tion without changing the value of the projected energy: 
The Lagrange parameter A is therefore superfluous. As 
a matter of fact, in a VAP approach one uses a Lagrange 
parameter only to speed up the convergence of the iter- 
ative procedure. 

It is important to realize that in a projected theory the 
only meaningful quantities are the projected ones. For 
example, the intrinsic density p{r) is not invariant un- 
der the transformations of Eq. This is simply the 
mathematical transcription of the fact that changing the 
particle number affects the intrinsic density. Conversely 
the projected density {f) is consistently invariant un- 
der the aforementioned transformation. 

In the demonstration above, Eqs. ([M)) - (|3ip we have as- 
sumed that the coefficients of the Bogoliubov transforma- 
tion (|12p are known. In the case of the full VAP approach, 
this is automatically the case because the U and V ma- 
trices are determined self-consistently by minimizing the 
projected energy, which is invariant under transforma- 
tions that change the number of particles. In the RVAP 
approach, however, to determine the U{S) and V{S) ma- 
trices one solves the standard HFB equations with a con- 
straint on the number of particles. The latter equations 
are obviously not invariant under the transformations of 
Eq. p4)) . This may generate a dependence of the RVAP 
solution on the Lagrange parameter A (or equivalently on 
(N)). Obviously 



d{m\H\m) 

d($((5)|iV|$((5)) 



-S. 



HFB 



(32) 

which illustrates that A„ can not be used to define the 
one neutron drip line in a PNP approach. It is interest- 
ing to realize that this dependence on A could be eventu- 
ally used to generate additional correlated wave functions 



I $(5, A)) in the RVAP M 
the projected energy [48 



aroach, thereby lowering further 
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FIG. 9: Upper panel: Projected energy for intrinsic wave 
functions with different average number of particles A'^ around 
the actual particle number A'^o = 182 in '^^^Os. This cor- 
responds to intrinsic wave functions with different Lagrange 
multipliers. Lower panel: corresponding Fermi level A. 

As an illustration we show in the upper panel of Fig. [5] 
the projected energy of ^^'^Os as a function of 6 for differ- 
ent values of the number of particles (or A) of the intrinsic 
wave function. As expected, we find that the minimum 
of the projected energy does not always correspond to 
the constraint (N) = A^Oi f^nd for a given constraint on 
the average particle number, the position of the mini- 
mum depends on S. In the lower panel the corresponding 
chemical potentials X{S) are plotted. 

If one restricts oneself to "one-dimensional" RVAP 
wave functions of the type |$((5)), it may happen, in par- 
ticular near the drip lines, that in the RVAP minimum, 
the underlying HFB wave function |<I'((5)} corresponds to 
a positive value of A. As emphasized before, this is with 
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no consequence since this A parameter does not define 
the drip hne. If one insists, however, in having a negative 
Fermi energy, it is always possible to slightly change the 
average number of particles of the intrinsic wave func- 
tion in such a way that A becomes negative, with the 
eventual cost of a small energy loss. In the illustrative 
case of ^^*0s displayed in Fig. [5] the energy cost is ap- 
proximately 15 keV to go from the RVAP minimum (at 
6 = 1.09) built on the HFB solution with average num- 
ber of particles (N) — No, to the RVAP minimum (at 
6 = 1.06) with (N) = No~ 1.0. The Fermi energy of 
the underlying HFB solution goes from +165 keV to -36 
keV. The drip line, defined from the two-neutron separa- 
tion energy S2n, remains unchanged. 



C. Nuclear halos and drip lines in a symmetry 
conserving approach 

In this section we investigate the effect of the parti- 
cle number projection on the size of the halo along the 
neutron drip line. As mentioned in the Introduction we 
should distinguish between halos in very light nuclei and 
in the heavier ones. We are aware that a mean field 
based approach may not contain enough correlations to 
describe the halo mechanism in very light nuclei. Nev- 
ertheless we shall first discuss the impact of the RVAP 
procedure on the archetypical case of halo nucleus ^^Li. 
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FIG. 10: Total projected energy (open squares) and neutron 
(plain circles) and mass (open triangles) r.m.s radius in ^^Li 
as a function of the RVAP parameter S. At each point 5 
the HFB solution is projected on good particle number. The 
minimum is attained ai S — 1.36. 

In the calculation of the nucleus "'^'^Li, with 3 pro- 
tons and 8 neutrons, the odd proton was treated in the 
equal- filling approximation (the lpi/2 state is the blocked 
state) and only the projection on neutron particle num- 
ber was carried out. Since the neutron number corre- 



sponds to a shell closure it is obvious that the HFB so- 
lution is not a super-fluid one. Figure fTOl shows the total 
projected energy and the neutron and mass r.m.s radius 
in ^^Li as function of the RVAP variational parameter 5. 
All calculations are done in the WS basis with the DIS 
interaction. The RVAP minimum always corresponds to 
a paired solution. In ^^Li, our original spherical HFB 
calculations with the DIS or Dl interactions do not pro- 
duce any halo. In fact, pairing correlations do not set in 
at all in this nucleus in the HFB calculations, even when 
the size of the box is increased up to 30 fm (thereby in- 
creasing the level density of continuum states). This is 
clearly viewed in Fig. [10] since the projected energy re- 
mains constant at = -47.48 MeV for 1.0<S < 1.18. 
In spite of multiplying the pairing field by the factor S 
during the iterations, pairing correlations are still identi- 
cally at convergence. Only for (5 > 1.18 do we observe 
the onset of significant pairing correlations. The total 
projected energy therefore decreases, continuum states 
begin to have a non-zero occupation probability, which 
contributes to the increase of the r.m.s neutron radius. 
At the minimum of the RVAP curve, both the neutron 
and mass radius have increased by about 2 %. The effect 
is marked but it is clearly not en oug h to reproduce the 
experimental halo in this nucleus [49l |. 
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FIG. 11: Neutron HFB density p^^^(r) in "Li at the HFB 
minimum (dotted line) and at the RVAP minimum (plain 
line). The dashed line shows the projected density p^^^(r) 
at the RVAP minimum. Calculations are done for the DIS 
interaction in the WS basis with Rbox = 20 fm. 

To better grasp the impact of particle-number projec- 
tion we show in Fig. [Tl]the neutron density in ^^Li in 3 
different cases. The dotted line corresponds to the stan- 
dard HFB calculation. The plain line corresponds to the 
density of the intrinsic HFB wave function |^'(i5)) at the 
RVAP minimum S — 1.36. We clearly see the formation 
of a " bump" which is a visual trademark of the nuclear 
halo. However, this solution is not physical since it is only 
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used to generate the variational space used in the RVAP 
procedure. Only the projected solution |$^(i5 = Smin)) 
at the minimum is physical. The corresponding projected 
density (dashed line) is slightly less extended than the 
underlying HFB solution. 



TABLE II: Table of two neutron drip line nuclei obtained 
using the DIS parametrization at the HFB and RVAP level. 
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2^°Rn 
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"«Ru 


136 Ru 


2^2Ra 


2^2Ra 


i«Pd 


140pj 


274^j^ 


274^j^ 


i52cd 


Cd 


278u 
284 p^ 
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Figures (TU] and [TT] suggest that the impact of particle- 
number projection may be instrumental in the forma- 
tion of sizeable halos, since the RVAP mechanism always 
guarantees a solution with non-zero pairing correlations. 
Since the chemical potential is irrelevant in a projected 
theory we should therefore, in principle, compute this 
quantity using projected energies and compare it from 
the results obtained using un-projected quantities. As 
emphasized earlier, the application of particle number 
projection in odd nuclei is not possible at the moment, 
hence the one neutron RVAP drip line is not accessible. 

We therefore carried out systematic RVAP calcula- 
tions of the two neutron separation energies, 5271) near 
the drip line using the DIS interaction. The particu- 
lar choice of the interaction is secondary in this study, 
since the focus is on the particular role of particle num- 
ber projection. The procedure was as follows: for a 
given drip line element {Z, N) from Table |T] the isotopes 
with — 4, N — 2, N and A^ -I- 2 neutrons were con- 
sidered. For each isotope, the RVAP procedure was 
carried with S = 1.0, 1.05, 1.50. The minimum of 
the RVAP curve was retained as the physical solution 
for every isotope. The two-neutron separation energy 
was calculated from the total RVAP-projected energies: 



S2n = BP^P{N,Z) - BP^P{N - 2,Z). The criterion 
S2n < was used to define the position of the new drip 
line. Table im shows the two-neutron drip line nuclei with 
and without the particle-number projection. These two 
drip lines differ by the isotopes of 6 elements: "^^Mgao, 
62Ca42, ^^Cr5o, i^2Sni22, 2«2Pti84 and ^""^P^ss- As we 
can read in the neutron number of these nuclei the differ- 
ences arise always close to the shell closures, where the 
pairing correlations are either very weak or vanishing. 

For all the elements located at the RVAP drip line, 
the Helm radius was computed, for the protons and the 
neutrons, based on the projected density p^^^{r). The 
quantity SR^o^'in) = R^i^oM - K^^in) obtained 
from these calculations is reported in Fig. [121 together 
with the original 5Rhaio{n) of the two un-projected drip 
lines (S'2n and Sn drip line). 
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FIG. 12: Measure of the halo: SRhaio{n) = 
Rgeomin) — RHeimin) foT RVAP-projected (plain squares), 
52n-unprojected (open circles) and Sn-unprojected (plain 
circles) drip lines. All results are based on spherical Gogny 
HFB calculations in the WS basis with the DIS interaction. 

The impact of particle-number projection is only sig- 
nificant in those nuclei that are unbound at the HFB 
level but bound in the RVAP-HFB. As can be seen from 
Table [m there are many nuclei that are particle-unstable 
{—Sin « A„ > 0) but two-particle-stable (5*2,1 < 0). In 
such cases, the halo is of course larger, sometimes signif- 
icantly larger like e.g. Cr or Fe, than the corresponding 
particle-stable isotope. Moreover the value of the halo 
calculated at the RVAP minimum closely follows the one 
calculated at the S2n drip line. The case of Cr is singular, 
in that the RVAP mechanism changes the two neutron 
drip line by 4 units, thereby considerably lowering the 
halo. 

However, beyond Z « 30, the differences between all 
approaches become relatively negligible. This goes along 
a very clear and definite trend towards smaller halos as 
the mass of the nucleus increases. Combining this obser- 
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vation with the fact that our mean-field approach, which 
includes the continuum and uses the best possible treat- 
ment of pairing correlations, fails to produce halos in 
light nuclei, it is tempting to conclude that halos are a 
trademark of few-body correlations only. As was recog- 
nized early, pairing correlations are a pre-requisite to the 
formation of halos in a mean-field approach indeed. How- 
ever, our work seems to further indicate that additional 
correlations beyond the symmetry conserving mean-field 
approximation are also mandatory. 

Figure [11] may suggest that for RVAP solutions the 
projected density profile is markedly different from the 
unprojected density. As mentioned already, only the pro- 
jected density bears a physical meaning by construction. 
The underlying HFB solution is only used to generate 
a set of highly pair-correlated projected wave functions. 
Nevertheless one may compare the behavior of SRhaio 
when using either the projected density p^^^{r) in the 
RVAP minimum or the underlying un-projected HFB 
density p^^^{r) in this same minimum to visualize the 
impact of projection itself. The difference in 5Rhaio is 
practically negligible (less than 0.01 fm) and certainly 
not on the same scale as differences coming from the in- 
teraction. 

We have also applied the RVAP formalism to the calcu- 
lation of neutron skins and we do not find any remarkable 
difference as compared with the HFB ones. This confirms 
the observation that neutron skins are, from a theoret- 
ical point of view, mostly sensitive to the details of the 
interaction (iso-scalar vs. iso- vector content), and from 
an experimental point of view, to the neutron excess but 
are not directly affected by the vicinity of the continuum. 

In conclusion, we applied our method to include contin- 
uum effects in spherical self-consistent HFB calculations 
with finite-range forces of the Gogny type to the case of 
nuclear halos and skins. Our calculations show that both 
the Dl and DIS parametrizations of the Gogny force lead 
to relatively small mean-field halos, that are of compara- 



ble size to most of the results obtained in Skyrme-HFB 
or relativistic Hartree-Bogoliubov theories. In particular, 
we do not find the giant halos in Neon and Zirconium 
isotopes that were reported in several publications. As a 
rule of thumb we observe that the size of the halo tends 
to decrease as the mass of the nucleus increases and only 
light nuclei feature decent-sized halos. By contrast, neu- 
tron skins are found to be very clearly related to the ratio 
(N-Z)/A. 

We also show that the impact of particle number pro- 
jection, before variation, is relatively important since it 
can change the position of the drip line. However, we find 
that particle-number projected continuum-coupled HFB 
theory, employing the most realistic form of the pairing 
interaction, cannot reproduce the large halos observed 
experimentally in very light nuclei such as ^^Li. This 
suggests a series of necessary conditions for a successful 
description of nuclear halos in the framework of mean- 
field theory: (i) the continuum must be properly included 
in the formalism (ii) the shell structure must be realis- 
tic enough (iii) pairing correlations must be present (iv) 
symmetry-breaking mean- field calculations, including all 
relevant deformation degrees of freedom, are probably 
mandatory (v) all such broken symmetries (in particular 
particle number) should then be restored (vi) probably 
configuration mixing such as e.g. GCM should also be 
included. 
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